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We present an analysis of diffusion in terms of the spontaneous density fluctuations in a non- 
tiiermal two-species fluid modeled by a lattice gas automaton. The power spectrum of the density 
correlation function is computed with statistical mechanical methods, analytically in the hydrody- 
namic limit, and numerically from a Boltzmann expression for shorter time and space scales. In 
particular we define an observable - the weighted difference of the species densities - whose fluctua- 
tion correlations yield the diffusive mode independently of the other modes so that the corresponding 
power spectrum provides a measure of diffusion dynamics solely. Automaton simulations are per- 
formed to obtain measurements of the spectral density over the complete range of wavelengths (from 
the microscopic scale to the macroscopic scale of the automaton universe). Comparison of the the- 
oretical results with the numerical experiments data yields the following results: (i) the spectral 
functions of the lattice gas fluctuations are in accordance with those of a classical 'non-thermal' 
fluid; (ii) the Landau-Placzek theory, obtained as the hydrodynamic limit of the Boltzmann theory, 
describes the spectra correctly in the long wavelength limit; (iii) at shorter wavelengths and at 
moderate densities the complete Boltzmann theory provides good agreement with the simulation 
data. These results offer convincing validation of lattice gas automata as a microscopic approach to 
diffusion phenomena in fluid systems. 

PACS numbers: Q5.20.Dd, 05.5Q.-fq, Q5.60.+W 

I. INTRODUCTION 

Frisch, Hasslacher and Pomeau pioneered a lattice gas automaton (FHP) as a microscopic model for incompressible 
fluids obeying the Navier-Stokes equation in two dimensions The FHP model was subsequently generalized to 
study diffusive phenomena in binary fluids using 'macroscopic' experiments [^-Q. Typically the observer would be 
interested in the evolution of the density proflle of 'red' particles in a system composed of 'red' and 'blue' particles 
where the color is a passive property used to distinguish species which otherwise do not differ from one another 
(it is necessary that they do in other circumstances Q). The diffusion coefficient is then evaluated by fitting the 
'experimental' profile to the solution of the diffusion equation subject to the appropriate boundary conditions . 

Because the FHP lattice gas lacks an independent coUisional invariant for energy, it is not suited for the modeling of 
thermal fluids. An appropriate generalization was realized by the construction of the model proposed by Grosflls, Boon 
^ ' and Lallemand (GBL) Their study was motivated by the analysis of the correlations of spontaneous fluctuations 
■ in lattice gas automata (LGA) in order to flnd whether the fluctuations power spectrum would be in accordance 
^ ' with those observed in actual fluids. Indeed the dynamical structure factor - the power spectrum of the density 
O fluctuations correlation function ~ gains its importance by providing insight to the dynamical behavior of the fluid 
[0, and the LGA was found to exhibit correct properties at global equilibrium: the spectra obtained by simulations 
of the GBL model present the same characteristics as those obtained from neutron- and light-scattering experiments 
in real fluids. In particular, in the hydrodynamic limit, one observes two shifted Brillouin peaks (corresponding to 
5^ ' the propagation of sound waves and their damping) and a central Rayleigh peak (corresponding to the diffusivity of 
\ entropy ffuctuations as a consequence of energy conservation). The GBL model was subsequently analyzed in detail 
by Grosflls et al^. 

The mixture of two real fluids exhibits a power spectrum in which the central peak is not a simple Lorentzian, 
even in the long wavelength limit 0. It has a spectral structure where it is diflicult to separate the contributions 
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from entropy fluctuations and from concentration fluctuations which are not decoupled in general (unless one of the 
two components is in trace amounts, in which case the two modes can be identified as they produce two independent 
central Lorentzians) . 

From the above considerations the idea emerged to analyze and measure the fluctuation correlations in a non- 
thermal two-species LGA fluid (in which the Rayleigh peak is absent) in order to study diffusion dynamics from a 
microscopic approach. In section || we present the model used for our numerical simulations. Section [II discusses the 
lattice Boltzmann theory for the analysis of the dynamical structure factor. The analytical results are developed in 
the hydrodynamic limit in section and are found to be in full agreement with the Landau-Placzek theory. In section 
^ we examine the different wavelength domains in terms of the Boltzmann propagator eigenvalues, and we present 
qualitative and quantitative analyzes of the theoretical results in comparison with the simulation data. The FHP 
spurious invariant is identifled and its effects are shown to be unimportant when the LGA is properly implemented. 
We conclude with some comments. 



II. THE MODEL 



The particles have unitary mass with no spatial extension and occupy the nodes of a triangular lattice with hexagonal 
symmetry. A particle can move along any of the six lattice directions (with unit velocity modulus) to one of the nearest 
sites or be at rest in its initial state (with zero velocity modulus). Particles interact via instantaneous local collisions 
which redistribute mass and momentum among the channels of each node at every time-step according to mass and 
momentum conservation. In a two-species system, particles are tagged either as 'red' or 'blue', and their color is 
redistributed randomly during the collisions independently of the mass redistribution; color is also conserved by the 
dynamics. The state of a node is given in terms of channel occupations. Here we use a description assigning a color 
to the channel. Since there are seven distinct velocities and two colors (one for each species), each node has seven 
pairs of channels. An exclusion principle is applied such that a pair of channels cannot be occupied by more than one 
particle (either red or blue) at any given time A corollary is that the equilibrium distribution takes the form of 
a Fermi-Dirac distribution ||l^ . The present formulation yields a convenient specification of the state of a node as a 
fourteen-bit word. 



III. THE BOLTZMANN FORMALISM 



The red mass density p^'^'^{r,t) is the number of red particles at node r at time t and the fluctuations 5p^'^'^{r,t) 
are defined in terms of the red channel occupations (i e {1, 6}, with 2b the total number of channels per node) 

b b 

5p-^''{v,t) = ^(5n,(r,i) = ^ [n,(r, i)- < n,(r, t) >] , (1) 

i=l i=l 

where <> denotes the equilibrium ensemble average; in basic equilibrium 

, , — r _ ( fS for i = 1, ... ,6 (red channels), , , 

< n,(r,t) >= - I j(i _ 5/) for i = 6 + 1, ... ,26 (blue channels), ^ ' 

with / the average density per pair of channels, and 6 the concentration of red particles. Note that 

26 

5]<n,(r,i) >= 6/0 + 6/(1 -0)-p-'' + p^'"^=p, (3) 

which defines the respective average densities per node. 

The 'red mass' dynamic structure factor 5''''^'^(k, w), defined as the space and time Fourier transform of the van 
Hove correlation function 

G'-^'^(r, t)^< 5p'-^'^{v, t) Sp^^'^iQ, 0) > , (4) 

is given by 
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where 5p'''^''(k, is the spatial Fourier transform of Sp'^^'^{r,t), and V is the total number of nodes, also interpreted 
as the volume of the lattice universe (here the lattice C is finite and has periodic boundary conditions). S''''^''(k, cu) is 
also expressed in terms of the kinetic propagator defined by 

T,j(k,t) = < 6n,(k,t)Sn*(k,0) > , = (7) 

where (Snj(k, t) is the spatial Fourier transform of Sni{r,t), and Kj = fj{l — fj). Using (|^), we write the dynamic 
structure factor (^) as 

oo b b 

predgred^^^^^^ ^ * ^ F,, (k, , (8) 

t— — oo i—1 j—1 

and the static structure factor (the Fourier transform of the equal-time van Hove function) as 

b b 

p-'^5-'^(k) = ^5]r.,(k,o)«, (9) 

1=1 j=i 

b b b 
i=l j=l j=l 



s'''"^{k) = 1 - fe . (11) 

We now evaluate the kinetic propagator in the Boltzmann approximation ||ll|. The lattice gas equation for the 
single particle distribution /i(r,t) reads |T2[ 

/,(r + c„t+l)-/,(r,t) + A({n,}). (12) 

Here A{{ni}) is the collision term, which is expanded around the stationary equilibrium distribution < Ui > to yield 

26 

A({< ni> + Sui}) = n.jSuj + 0{SnjSnk), (13) 

where we have used the property A({< Ui >}) = which follows from mass conservation. The explicit form of Qij is 
given in terms of the transition matrix A{s s') between pre- and post-coUisional states s and s' respectively 

- ^ A{s-.s){s^- s,)s, <„ ■ >< 1 - „ > ■ (14) 

This result is obtained with the assumption that particles on different channels of the same node arc uncorrelated 
before collision (Boltzmann ansatz), i.e. by factorizing the averages < Uj > {i ^ j). Combining (|l^) and (p^, we 
obtain the linearized Boltzmann equation which reads in k-space 

2b 

6n,{k,t+l) = ^e-'''"^' {6,,+n,,)dnj{k,t). (15) 

i=i 

Equation ( p^ is straightforwardly solved by iteration; inserting its solution into (^) yields 

r,,(k, t) = [e-'''"^ ■ (5 + fl)]* . , (t > 0). (16) 
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Here [e-'^-'^y = fijie'^^"^ 



is a diagonal matrix. From (||) and (|T^) we obtain 



b b 

i=i j=i 



1 



1 



(17) 
(18) 



where TZe denotes the real part. This expression for the dynamic structure factor is exact within the Boltzmann 
approximation, but the explicit analytical inversion of the 6x6 matrix in (|l^) cannot be performed in all generality. 
However perturbation methods can be used to compute analytically S^'^'^(k,uj) in the hydrodynamic limit: |k| — > 
and uj ^ 0(|k|), Odkp) (Section IV). Beyond the long wavelength - long time domain, one has recourse to numerical 
evaluation of ( p^ ) to compute the 'Boltzmann' power spectrum (Section ^). 



IV. THE HYDRODYNAMIC LIMIT 

A. The hydrodynamic modes 

We first notice that the linearized collision operator 17 is not symmetrical, with the consequence that its left and 
right eigenvectors are not each others transpose. However when the detailed balance condition is satisfied, the matrix 
product ilijKj is symmetrical ||l^, and the left and right eigenvectors of fl are related by \(f))i = Ki it can also 
be shown that to each of the N collisional invariants corresponds an eigenvector (A„|(n = 1, A^) of fl, with zero 
eigenvalue. The components are given by the conserved quantities carried by each channel i : 
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(19) 



From these considerations and by analogy with the thermal scalar product introduced in |8|, we define the colored 
scalar product 

6 

{A\B) =Y,Mc^)l^^B{c,), (20) 

i=l 

where the weight Ki depends on density and concentration. S symmetrical matrix, the colored scalar 

product has the symmetry 

26 2b 

{A\n\B) = {B\n\A) = ^^A(c,)f^.,A^,i?(c,). (21) 

i=l j=l 

Following a method introduced by Resibois [Q, we consider, as the starting point, the propagator ( pj| ) which is 
the t-th power of the non-symmetrical matrix e^^ *^ • ((5 + r2), and we use the eigenvalue problem formulation 

e-*- -{d + n) IV'p(k)) = e^-^'') IV'p(k)) , (22) 
(0^(k)| e-*- ■id + n)= e^-^'') (0^(k)| . (23) 

The eigenmodes of the propagator may be separated into two groups: the slow modes, corresponding to eigenvalues 
Zfj,{k) close to zero when fc(= \k\) tends to zero, and the fast modes corresponding to eigenvalues TZe Zfj,(k) < leading 
to exponentially fast decay. The latter are the kinetic modes; the slow modes which decay infinitely slowly when 
fc — > will be identified as the hydrodynamic modes. They are the dominant modes in the hydrodynamic regime 
where the kinetic modes can be neglected. 

For |k| = the matrix e"'*' '^ • (5 + O) reduces to S + whose eigenspace spanned by the eigenvectors ( |l9| ) has the 
dimension given by the number of collisional invariants (here 4). For |k| 7^ but small, we can express the eigenvectors 
of e"*'' '^ • (5 + $7) as a linear combination of the collisional invariants, and we can expand e"*'^ '^, |'0^(k)), (0^(k)|, 
and z^(k) respectively as 
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(24) 



with • 



= (5,;, c, i,j — 1, 26, where ^ denotes the projection of onto k. Substitution of the first and second 

(25) 
(26) 



expressions of ( |24| ) into (|2^) and identification of the successive powers of k yields the hierarchy 
Oik'): n|^W) = (c,+z«j)|^f ), 



Oik'): n|^f) = (c,+zW 5)141)) + 
The solution to zero-th order is straightforward; one has 

N 



(27) 



(28) 



where the coefficients 6„ are to be determined subsequently. 

The first order solution is obtained by taking the scalar product of {A^ \ with ( p6| ) where the previous order solution 
is substituted; the result has the form as an A^-dimensional eigenvalue problem: 



AT 

E(A„| (c, +z«5) |A„)6„, =0, 



(29) 



which yields the four eigenvectors \ipli'^) and eigenvalues 



.(1) 



shear mode: 
acoustic modes 



IV'f ) = \p±) 

IVl")) = \P,)+Cs\M) 
color diffusion mode: I'tpim-) — i^b\R) ^ Kr\B) 



0, 

-Cs , 

0. 



(30) 



Here \M) is the sum of \R) and and P± are the projections of the momentum onto k and perpendicular to k 

respectively, Kj = - fj) with j = 1, for Kr, and j = 6 + 1, ...,26 for Kb, and Cs = i(P,\P,) / {M\M))'/' will be 
identified as the speed of sound (here Cs — ^3/7). 
We define the currents \jfj,) as 



(31) 



and we note, by multiplication of ( p6| ) by {Tp^\^'^\ that the currents are orthogonal to {^p^J'^^l- As a consequence the 
currents do not belong to the f2 kernel, and we may write the formal solution to the first order equation (Eq) as 



1 ^ 



(32) 



The coefficients 6^^ are determined by substituting IV'^i^'') by ( |32D in the second order equation (|27|), and multiplying 
the result by {■)pi°'' \ iv ^ fi) to obtain 



6,.(4°)|4°)}(zW-z, 



The expression for zj^^ follows from the evaluation of the product of equation (^) with {ipjj^^ \ , which yields 



= 



(33) 



(34) 
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We anticipate that zf^ is the kinematic viscosity (j/), that zf = z(2) is the sound damping (F), and that z^^X is the 
color-difFusivity (-D), as wih be justified subsequently by the analysis of the power spectrum. 

Explicit evaluation of (33) shows that the only non-zero off-diagonal elements of the matrix formed by the h^iv^ are 
the two coefficients 



= - 6_ , = 



2 c, 



(35) 



The diagonal elements remain unknown, but this is unimportant because, as will be seen, they do not contribute 
to the power spectrum (Section ^) . 

We have now identified the four hydrodynamic modes in the LGA. The shear mode and the acoustic modes are 
independent of color related properties, and the mode |i/'dift) describes color diffusion only. As will be shown below, 
the density power spectrum reflects this property. Notice that the purely diffusive behavior of color is related to an 
observable (defined below) which is neither the concentration of one of the components nor the difference between 
the two concentrations HH. 



B. The dynamic structure factor 



Transposing Eq.(E3h which defines the left eigenvectors of e • (^ + rj), and multiplying the result by e on 



the left, wc find that ((/)^| and j?/;^) are related by 



m,, 



(36) 



where is a normalization constant. If the eigenvectors and |'(/'^) form a complete bi-orthonormal set, i.e. 



26 



^ |^^(k)}(</)^(k)| = 6 , and ((/.^(k)|^,(k)) = 8^, , 



(37) 



we may write e • (^ -I- f2) as 



26 



e-'*-- • (,5 + f2) = ^ |^^(k))e^-W(0^(k)| 



(38) 



We will use this expression to recast the spectral function i^'"'^''(k, w) ; we first rewrite ( |18D as 

6 6 r 



i=\ j=i 

{R\ 



1 



1 1 



= {R\ 



1 



S\R) 



{R\ i + ld\R) , 



{s + ny^ + l-S\R) 



(39) 



where the last equality is obtained by noticing that the action of (S + ft) upon \R) is the identity operation since 
\R) belongs to the kernel of O. Then making use of (38), we find 



F^^\k,cu) = J2{R\^,{k)) 



1 



giw-Zj,(k) _ I 2 

with f^^^ {R\ij^{k)){(j)^{k)\R) , and = 



,{k)\R) = 5]A/;2?^, 
1 1 



(40) 



We observe that each mode fi contributes a spectral line whose amplitude depends on 2?^;. This factor becomes large 
for {iu — z^) 0, that is for small z^ in the limit of small oj. The modes for which z^(k) tends to zero at long 
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wavelength are precisely the slow modes identified in (|30|). So we may approximate (^) by neglecting the fast kinetic 
modes in the sum over /i. It is then consistent to make use of the approximation (e^ — + i « + 0{x) for 
a; ^ 1 in the evaluation of V^; with (]2^), we obtain 

+ (41) 



The final step is the evaluation of Af^ in terms of the /c-expansion of and \tpn)] this is accomplished by expressing 

i: 

orders 0{k'^) and O(k^) respectively, we find 



{f"*! in terms of \ipl?''), using ( ^41) in (|36|), expanding in powers of k, and identifying the successive orders. To 



(43) 

{''Pi, Wii ) 



which results are inserted into (|o|) to yield 



We now discuss the evaluation of for each hydrodynamic mode separately 

^1° 



• /i =± . As the vectors {R\ and are orthogonal, 7V_l = 0, and the shear mode will not show up in the density 
fluctuations power spectrum. 

• /i = ±. To order 0{k^), Af± ~ ^^^/{nr + '*&)• The computation of the next order requires, in principle, complete 
knowledge of l^/'^^''), but in fact this is unnecessary because the terms including the unknown fe^^ cancel each other; 
thus we obtain A/± = + i^b)[^ T ikT/{2cs) + ...). 

• /i = diff. To order 0(/c°), TVdiff — bKrKii/{Kr + Kb). The contributions to order 0{k^) all vanish identically; so the 
expression to 0{k^) is correct up to corrections of 0{k^). 

Combining Af^ with the explicit expression of Vf, (using (1^), (|l|) and ^ in (|l])), we obtain F''^'^ (k, lu) , and 
therefrom the analytical expression of the dynamic structure factor which reads 



^ ' ^ K-., -I- ^ 



Tk Csk±Lu 



V(^±csfc)2 + (rP)2 c, (w±c.fc)2 + (rF)^ 

2DP 



(46) 



At fixed value of k, the spectrum consists of a Brillouin-doublet centered around ifcCg, and of a central peak charac- 
terizing color diffusion. To first order in fc, Xm{z^) yields the frequency shift of the spectral peak corresponding to the 
propagating modes (/i = ±), and to second order in fc, 7?.e(z^(fc)) yields the dissipation coefficients which determine 
the hne-widths of the spectral components. Note that (|^) and (0) yield J^^ dujS'''"^{k,Lu)/27rS'''"^{k) = 1; the 
spectra shown in the figures are normalized accordingly. 
Considering the fluctuations 

5pd.ff(r, t) = p<j.tt(r, t)- < pd.ff(r, t) > (47) 

of the observable pditt deflned as 

— „red '^r blue //|n\ 
Pditt — P P , (,4»J 
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the corresponding spectral density 

oo 

reC t=~oo 

can be computed straightforwardly along the lines of the evaluation of the power-spectrum 5^'^'^(k, cj); since I'tpi^l-) is 
orthogonal to the other eigenvectors, there is no coupling with the other modes, and the dynamic structure factor is 
a single Lorentzian 

with similar normalization as for S^'"^(k,uj). The power spectrum ( pO| ) characterizes color diffusion alone, a feature 
which will be used in the analysis of the simulation data in the next section. 



V. THE POWER SPECTRUM 



The above results obtained in the hydrodynamic limit are in accordance with the Landau-Placzek theory for 
continuous fluids |^]. Now the hydrodynamic theory breaks down at short wavelengths, but the Boltzmann theory 
should remain valid down to k values where the kinetic domain starts. Lattice gas automata are appropriate model 
systems to investigate quantitatively the various regimes covering a wide range of wavelengths. In order to characterize 
the wavelength domain of the various regimes, we consider the quantities k£f, where is the mean free path, 

and / — p/pmax, the reduced density (or the average density per channel). Accordingly the hydrodynamic regime is 
defined by k£f ^ 1, the generalized hydrodynamic regime (Boltzmann regime) hy k£f < 1, and the kinetic regime by 
kif > 1. 

We now discuss the power spectra obtained from automaton simulation data and we compare the results to the 
analytical Landau-Placzek expressions and to the predictions of the lattice Boltzmann theory. For the latter, we use 
the eigenvalue spectrum of the propagator F which can be evaluated numerically over the complete fc-domain, so 
extending the computation of the power spectrum to the region of k values where analytical evaluation can no longer 
be performed. In Fig.l we show a typical eigenvalue spectrum as computed numerically, where the fourteen modes of 
the 14-bit model described in Sec.^ can be distinguished. 

The numerical experiments are performed using a 'color' FIIP-3 lattice gas (see Sec.||) at equilibrium. The lattice 
size is 256 x 256 nodes, and the simulation duration is 40.000 time steps. Spatial Fourier transforms are computed 
at every time-step, and timc-Fourier transforms are taken over intervals of 16384 time-steps shifted by 20 time-steps 
for averaging; data shown in the figures are smoothed by low-pass frequency filtering. 

In Fig.l. a we observe that in the range < fc < 0.4, corresponding to wavelengths A > 15^*0 (with Iq the lattice 
unit length), the four slow modes are well separated from the kinetic modes, and their behavior is correctly given by 
the hydrodynamic expressions. We therefore expect that in this domain, the Landau-Placzek theory should provide 
a correct description of the dynamic structure factor. This is indeed confirmed by the comparison between the 
results obtained from simulation data and the theoretical predictions as shown in Fig. 2, where we also notice that the 
Landau-Placzek spectrum is indistinguishable from the Boltzmann spectrum. 

The diffusion coefficient D{f, 9) is obtained straightforwardly from the diffusive mode power spectrum: S'^''^f{'k, cj) 
is a single Lorentzian ( po| ) with half- width Aui — Dk^; so plotting Auj as a function of fc^ (see Fig. 3. a) yields, by 
least-squares fit, a slope whose value provides an experimental measure of D. In Figure 3.b, we show the diffusion 
coefficient as a function of density: we observe that the exact Boltzmann result [|l^ is in good agreement with the 
lattice gas simulation data up to / « 0.25; for larger / the theoretical prediction deviates progressively from the 
measured values indicating that the molecular chaos assumption becomes invalid at high densities. 

As k increases from 0.4 to 1.4, there is still a distinct scale separation between slow and fast modes (see Fig.l), but 
the eigenvalues of the slow modes depart significantly from the hydrodynamic prediction, indicating the breakdown of 
the local response hypothesis: the transport coefficients become fc-dependent. As a result, the Landau-Placzek theory 
does no longer describe the power spectrum correctly - for instance, there is a noticable spectral line broadening in 
S'^-'^f f [k,uj) (see Fig. 4. a) -, but the complete Boltzmann spectrum is in good agreement with the simulation results, 
as seen in Figs. 4. a and 4.b. We have also observed that even at rather short wavelength (A ^ lO^o) the experimental 
data can still be approximated with a Landau-Placzek spectral function if the transport coefficients are parameterized, 
showing that a hydrodynamic type description holds qualitatively down to quite short wavelengths. For k > 1.5, all 
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modes exhibit comparable decay rates (see Fig.l.a), and all modes with even parity in c± contribute significantly to 
the power spectrum. In this domain, the Landau-Placzek theory is invalid and the Boltzmann computation provides 
good agreement with the experimental spectrum (down to A « Mq) as shown in Fig.4.b. 

From the agreement between the experimental data and the Boltzmann results, we found that the Boltzmann theory 
remains valid up to reduced densities oi f ~ 0.25. At higher densities the discrepancy between the Boltzmann spectral 
density and the experimental power spectrum (see Fig. 5. a) reflects the failure of the Boltzmann theory to evaluate 
correctly the transport coefficients. Here the contributions of 'ring-collisions' should be included in the evaluation 
of the diffusion coefficient. Finally we note that at very high densities (/ = .9), we observe a slight coupling between 
the color diffusion mode and the sound propagation modes as illustrated in Fig.5.b; so far we have no theoretical 
analysis of this effect. 

We close this section with some comments on the spurious invariants which have been known to plague lattice gas 
automata with spurious conservation laws The present model is not free of this 'contamination': figure 6. a shows 
clearly that TZe{z±{k)) (i.e. le^-L^'^)] ^ 1) as fc ^ 27r (which corresponds to a wavelength of one lattice-unit), 
denoting a mode which persists once it is excited. The corresponding spurious invariance can be interpreted as the 
conservation of total transverse momentum on even and odd lines of the lattice every two time-steps {Xm{z±{k)) tt, 
and therefore e^-^'^'^-' ^ — 1 ). Consequently it is important to choose initial conditions such that the total transverse 
momentum on odd and even lines be both rigorously zero. This may easily be realized by implementing the particles 
pairwise with opposite velocities on the same node. Without this precaution, the power spectrum may exhibit spurious 
effects as illustrated in Fig.6.b. The full line shows the correct spectrum. Now if we denote the reference direction by 
x, and measure the spectrum with k oriented along the y-direction (orthogonal to the x-axis), we obtain the spectrum 
represented by the dotted line in Fig.6.b: the speed of sound is not affected, but the line-width of the shifted peaks is 
incorrect, and the diffusive mode is propagative, despite the fact that the spectrum is measured in the long wave-length 
domain ( A « 64^o). With k oriented at an angle of 30 deg. off the a;-axis, we observe (see dot-dashed spectrum in 
Fig.6.b) that there is a significant shift in the speed of sound, that the diffusive mode is practically absent, and that 
the kinetic modes which couple to the transverse momentum invade the spectrum. 

VI. COMMENTS 

We have presented a two species non-thermal lattice gas automaton for which we have developed the lattice Boltz- 
mann theory and performed automaton simulations. The emphasis is on the fluctuation correlations in order to 
obtain a microscopic analysis of diffusion dynamics in contrast to earlier studies based on macroscopic approaches. 
On large space- and time-scales we find spectral features of the dynamic structure factor in accordance with those of 
real fluids described by the Landau-Placzek theory. Because of the intrinsic simplicity of the lattice gas model, the 
various wavelength regimes can be easily identified, and the propagator spectrum can be used to to compute the power 
spectrum over the full wavenumber domain, and to test the the validity of the Boltzmann hypothesis. The present 
study based on the analysis of spontaneous fluctuations offers a microscopic approach to diffusion , and, through the 
identification of a purely diffusive mode associated to color transport, complements and supports earlier macroscopic 
investigations. 
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FIG. 1. Eigenvalue spectrum of the 14-bit model propagator: Boltzmann computation (full lines) and hydrodynamic limit 
(dashed lines). The reduced density and concentration are / = 0.15 and 9red = 30% respectively. The wave-number k — |k| is 
given in reciprocal lattice units 



FIG. 2. Spectral measurement of the diffusion coefficient: (a) Alo — Dk^ , f = 0.3 and 9red = 30% (open squares), 50% (black 
dots); the least-squares fits (solid line) coincide, (b) D — D{f,9red), simulation data (black dots) and Boltzmann prediction 
(solid curve); the size of the black dots corresponds to the largest error bar {\AD/D\ < 2%). 



FIG. 3. Power spectra of 'red' density fluctuations (a) and of pdiff fluctuations (b) at low density and small k. Comparison 
of experimental data (full lines) with theoretical predictions: the Boltzmann results and the Landau-Placzek spectra coincide 
(dashed lines). Density / = 0.15; concentration 9red ~ 30%; wave number |k| = 0.098 reciprocal lattice units; u is given in 
reciprocal time units {2n/T, where T is total number of time-steps); the spectral functions are given in reciprocal u units. 



FIG. 4. Power spectra at low density and high k: |k| = 0.49 (a), i.e. A « 12^o; Ik[ — 1.57 (b), i.e. A « 4£o. Experimental 
data (full lines), Boltzmann spectrum (dashed lines), Landau-Placzek theory (doted lines). / = 0.15; 9red = 30%. 



FIG. 5. Spectra at high density: / — 0.5 (a) and / — 0.9 (b). Experimental data (full lines) and Boltzmann prediction 
(dashed line). Concentration 9red ~ 30%; wave number |k| = 0.098. 



FIG. 6. Effect of spurious invariant, (a) Boltzmann propagator eigenvalues (full lines) with spurious invariant (dashed line); 
(b) shows the results of two simulations: the spectrum plotted as a full line is obtained with initial conditions for which the 
oscillating transverse mode is not excited; the two other power spectra are obtained for two k's with different orientations using 
'incorrect' initial conditions (see text). Density / — 0.10; concentration 9red ~ 50%; wave number |k| = 0.098. 
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